function [ v ] = asymptoticvbar( sigma, N )

x = load('tableData.dat');

H = x(:,2) ./ 2;
N = x(:,1);

omega = 4*pi/3;

epsilon = sqrt(4*sigma./N);


v = omega./(4.*epsilon.*N) .* ( 1 + epsilon/pi .* log(2./epsilon) + epsilon./pi.*(-9.*N/5 + 2.*(N-2).*log(2) + 3/2 + 4./N .* H));

end

